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ABSTRACT 


The  energy  gap  of  Sn-Sn  and  Sn-Pb  tunnel  junctions  has  been  measured 
with  accuracy  better  than  one  per  cent.  The  values  of  the  energy  gaps 
at  0°K  are  found  to  be  2Agn(0)  =  1.22  mv  and  2Ap^(0)  =  2.6  mv  .  The 
ratio  2Agn(0)/KgT£  -  3.68.  It  is  experimentally  confirmed  that  the 
variation  of  [A(T)/A(0)]2  is  a  linear  function  of  t  =  T/Tc  just  below 
Tc,  as  indicated  by  the  BCS  theory.  The  transition  temperature  of  Sn 
films  is  (3.83  +  0.01)°K.  The  relative  jump  in  the  specific  heat  at 
Tc  is  calculated  from  the  I  -  v  traces  to  be  (1.61  +  0.02). 

The  phonon  structure  in  Sn  and  Pb  has  been  investigated  at  2°K 
and  4.2°K  respectively. 
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CHAPTER  1 


INTRODUCTION 

One  of  the  most  important  facts  about  the  nature  of  superconductivity 
is  the  existence  of  a  gap  in  the  energy  spectrum  of  a  superconductor. 

Many  techniques  are  now  available  for  measuring  such  a  gap.  Electron 
tunneling,  discovered  by  Giaever  (1961),  is  one  of  the  most  direct 
techniques  to  investigate  many  of  the  superconducting  properties.  This 
technique  has  been  initiated  in  this  laboratory  by  Adler  (1962)  and 
Rogers  (1963),  then  extended  by  Penner,  Khanna,  and  Keeler.  The  ultra¬ 
sonic  attenuation  experiment  is  also  a  very  useful  tool  to  investigate 
the  energy  gap  anisotropy  and  its  temperature  dependence. 

This  thesis  contains  a  report  of  the  first  determination  of  the 
relative  jump  in  electronic  specific  heat  of  Sn  at  the  transition  tempera¬ 
ture  Tc,  from  the  current-voltage  traces.  This  measurement  might  be 
useful  especially  for  elements  for  which  the  electronic  part  is  a  small 
fraction  of  the  total  specific  heat.  The  energy  gap  of  Sn  as  a  function 
of  temperature,  calculated  from  ultrasonic  attenuation  data  and  directly 
measured  in  an  electron  tunneling  experiment,  is  compared  with  the  BCS 
theory. 

The  electron  tunneling  is  not  a  good  method  of  exploring  the  aniso¬ 
tropic  properties  of  the  superconductor.  With  specimens  composed  of 
polycrystalline  films,  one  can  only  obtain  the  average  of  a  certain 
measurable  quantity  over  the  area  of  the  tunnel  junction. 
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CHAPTER  2 


MICROSCOPIC  THEORY  OF  SUPERCONDUCTIVITY 
A.  The  BCS  Theory 

Kamerlingh  Onnes  (1911)  investigated  the  temperature  dependence  of 
the  electrical  resistivity  of  a  mercury  sample  at  low  temperatures.  He 
found  that  at  about  4°K,  the  resistance  falls  abruptly  to  zero.  This 
extraordinary  phenomena  he  called  superconductivity,  and  the  temperature 
at  which  it  appears,  the  critical  temperature  Tc.  Below  Tc  the  super¬ 
conducting  behaviour  can  be  quenched  and  normal  conductivity  restored 
by  the  application  of  an  external  magnetic  field.  It  was  not  until 
(1933) ,  that  Meissner  and  Ochsenfeld  investigated  the  magnetic  properties 
of  superconductors  in  some  detail.  They  discovered  that  the  magnetic 
induction  in  a  superconductor  vanishes  (Meissner  effect) .  The  discovery 
of  the  Meissner  effect  was  required  to  establish  that  the  superconductive 
state  is  a  true  thermodynamic  phase.  In  view  of  the  Meissner  effect, 

Gorter  and  Casimir  (1934)  were  able  to  develop  a  full  treatment  of  the 
superconducting  phase  transition.  Their  treatment  indicates  that  the 
transition  is  of  second  order.  Consequently,  there  is  a  jump  in  the 
specific  heat  at  T  =  Tc.  The  magnitude  of  the  jump  is  twice  the  electronic 
specific  heat  of  the  normal  metal  at  the  same  temperature.  Also  below  Tc 
the  entropy  of  the  superconducting  phase  is  lower  than  that  of  the  normal 
one,  indicating  that  the  former  is  the  state  of  higher  order.  This  ordering 
will  later  be  shown  to  follow  from  a  condensation  of  electrons  in  momentum 
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space. 

Application  of  a  magnetic  field  greater  than  Ec  the  critical  field 
permits  a  determination  of  the  energy  difference  between  the  normal  and 
superconducting  states.  This  difference,  known  as  the  condensation  energy 
is  extremely  small,  of  the  order  of  10-8  ev  per  atom.  Among  other  diffi¬ 
culties,  this  small  energy  difference  stands  as  a  basic  one  in  constructing 
a  microscopic  theory  for  superconductivity. 

A  first  successful  theory  has  been  developed  by  Bardeen,  Cooper,  and 
Schrieffer  (1957);  (known  as  BCS  theory).  The  basis  of  the  theory  was 
provided  by  Cooper  (1956)  who  showed  that  in  the  presence  of  a  Fermi  sea 
of  electrons  if  two  electrons  had  an  energy  slightly  greater  than  the  Fermi 
energy  then  in  the  presence  of  an  attractive  interaction  they  would  form  a 
bound  "pair”  state.  BCS  showed  that  in  the  presence  of  an  attractive 
interaction,  the  electrons  in  the  neighbourhood  of  the  Fermi  surface 
condense  into  a  state  of  lower  energy  in  which  each  electron  is  paired 
with  one  of  opposite  momentum  and  spin.  The  theory  predicts  the  existence 
of  an  energy  gap  which  separates  the  state  with  the  largest  possible  number 
of  Cooper  pairs  from  the  state  with  one  pair  less.  This  leads  to  the 
correct  thermal  and  electromagnetic  properties  to  display  superconductivity. 

The  attraction  between  electrons  necessary  to  form  Cooper  pairs  can 
in  principle  be  due  to  any  suitable  kind  of  interaction.  The  discovery 
by  Maxwell  (1950),  that  for  many  superconductors  the  critical  temperature 
Tc  depends  on  the  isotopic  mass,  showed  that  for  these  substances  the 
attractive  interaction  is  one  between  the  electrons  and  the  lattice.  The 
simplest  way  of  calculating  this  attractive  interaction  is  to  describe  it 
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as  the  emission  of  a  virtual  phonon  by  one  electron,  and  its  absorption 
by  the  other.  Consider  the  processes  defined  in  fig.  (2-1),  in  which 
(a)  an  electron  in  a  state  k  emits  a  phonon  q,  and  is  scattered  into 
state  k  -  K;  the  electron  in  state  k  absorbs  this  phonon  and  is  scattered 


into  k'  +  K.  (b)  A  phonon  -q  is  emitted  by  K''  ,  before  being  absorbed  by 
k.  The  matrix  element  for  the  transition  corresponding  to  these  processes 
is  given  by 


(2-1) 


<k  -  K,  k  +  K|V 


->  -k 

k,  k 


>  = 


k-K1 

(e(£)  -  e(li-K)}2- («»>_,.) 2 


In  this  expression  ^  ^  is  the  matrix  element  for  electron-phonon 
interaction,  e(k)  is  the  Bloch  energy  of  an  electron  in  a  state  k 
measured  relative  to  the  Fermi  surface,  and  iTu)_>.  is  the  phonon  energy, 

q 

which  is  of  the  same  order  as  the  Debye  energy  KB0B.  It  is  obvious 
that  the  interaction  is  attractive  (negative)  only  for  the  region  of 
energies  where 

(2-2)  |e(£)  -  e  (k-K)  |  <  -ha)^  ~  kB0D 

q 


and  repulsive  otherwise.  In  this  range,  the  matrix  element  for  the 
interaction  can  be  written  as 


(2-3) 


-  2 


!**£.  k- 

•fio) 


K1 


q 


The  corresponding  matrix  element  for  the  screened  Coulomb  repulsion  between 
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Fig.  (2-1) 


Schematic  representation  of  electron-electron  interaction 


by  the  exchange  of  phonons. 
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electrons  would  be 


(2-4) 


4ne2 

K2  +  X2 


where  X  is  the  screening  parameter.  The  criterion  for  superconductivity 
to  occur  is  for 


(2-5) 


-  V 


M+ 


<  -  2 


k,  k-K1 


q 


4lle 
K2  +  X2 


to  be  negative.  The  average  is  taken  over  electron  states  which  satisfy 
the  inequality  (2-2). 

Because  of  the  Fermi  statistics  of  the  electrons,  the  matrix  element 
for  all  possible  interactions  which  take  a  pair  of  electrons  from  any  two 
k-values  to  any  two  others,  alternate  in  sign  and  they  are  all  of  roughly 
equal  magnitude.  The  result  is  a  total  interaction  energy  that  is  smaller 
than  the  individual  terms.  In  order  to  have  matrix  elements  of  a  single 
sign,  BCS  had  to  associate  all  possible  k  values  in  pairs  (k,  k  ),  and 
require  that  either  both  or  neither  member  of  a  pair  be  occupied.  The 
conservation  of  momentum  requires  that  k  +  k  =  P  should  be  the  same 
for  all  bound  pairs.  However,  it  is  energetically  most  favourable  to 
restrict  the  pairs  to  those  for  which  P  =  0  and  spin  is  zero.  BCS 
proceed  to  calculate  the  superconducting  ground  state  energy  entirely 
from  the  correlation  between  paired  states  (kt ,  -kl),  by  phonon  and  screened 
Coulomb  interactions.  A  pair  will  be  designated  by  a  wave  vector  k 
independent  of  spin. 
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The  Hamiltonian  for  the  ground  state  can  be  expressed  as 


(2-6)  Hred  =2  l 


-> 

k 


e  b*  b 
k  k  k 


l 


k,k 


b*  b 
kk  k  k 


where 


C*  € 
kf  -ki 


k 


C  C 
-kl  kt 


are  the  creation  and  annihilation  operators  for  pairs,  C^_  and 

ka  ka 

are  the  single  electron  state  creation  and  annihilation  operators, 

V  is  the  strength  of  the  attractive  interaction  of  one  pair.  The 

S' 

introduction  of  pair  operators  removes  all  the  singly  occupied  states 


from  the  Hamiltonian  Hrecj. 

The  ground  state  wavef unction  of  Hre(^  is  given  by 


(2-7) 


I  \b>  =  II  (u  +  v  b  )  nb  > 

k  k  k  k  u 


where  v  2  is  the  probability  that  the  k  ——  pair  is  occupied,  u  2  is 
k  k 

the  probability  that  it  is  not  occupied,  and  |  tJ)  >  is  vacuum  state, 

The  normalization  condition  requires  that 


2  ,  2 

u  +  v 

k  k 


=  1  . 


The  condensation  energy  corresponding  to  Hred  for  the  state  (2-7)  is 


(2-8) 


wc  =  2  I 


£  V 
->  -> 

t  k  k 
k 


I 


kk 


v  u  v  u  V 

->  ->  >.  ->  -» 
k  k  k  k  kk" 


The  first  term  gives  the  difference  of  the  kinetic  energy  between  the 
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superconducting  and  the  normal  phases;  the  factor  2  arises  because  of 
the  pairing  states.  This  first  tern  is  smaller  than  the  second  term, 
which  gives  the  correlation  energy  for  all  possible  transitions  from  a 


pair  state  k  to  another  k  .  For  such  a  transition  to  be  possible,  k 


7^"  * 


must  initially  be  occupied  and  k  must  be  empty.  The  simultaneous 


2  2 

probability  of  this  is  given  by  v  u  .  The  final  state  must  have 

k  %■'  2  2 

k  empty  and  k  occupied;  this  has  probability  v  u^  .  The  total 

k  k 

correlation  energy  is  given  by  the  sum  of  the  product  of  the  matrix 

element  V  by  the  amplitudes  v  u  and  v  u  over  all  the 
kk  ^  k  k  k  k 

possible  values  of  k  and  k  . 

The  necessary  condition  for  W  to  be  minimum  is 


(2-9) 


v 


-> 

k 


(1  - 


£-> 

__k 

E+ 

k 


(2-10) 


u 


S 


1 

i  (1  +  Js ) 

z  E 

-> 

k 


where 


(2-11) 


-> 

k 


=  (e 


k 


+ 


and 


(2-12) 


-> 

k 


-  I 


V  A 

-  y—y  *  ~y  * 

kk  k 


S' 


2  E 


S' 


Figs.  (2.2a)  and  (2. 2b)  illustrate  equations (2-9)  and  (2-11)  respectively. 

BCS  make  use  of  the  simplifying  assump tions  that  the  Fermi  surface 
is  a  spherical  one,  and  V  is  constant  for  electrons  m  the  narrow 

SS" 
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region  defined  by  (2-2)  and  vanishes  elsewhere.  Assuming  symmetry  of 
states  on  either  side  of  the  Fermi  surface  (ep  =  0),  and  introducing 
N(0)  the  density  of  single  electron  states  of  one  spin  in  the  normal 
state  at  £  =  0,  equation  (2-12)  becomes 


(2-13) 


N(0)  V 


■fiOLy 

=  /  q 

0 


-  KBeD 


de 


42  +  A2(0) 


The  limit  of  integration  is  the  phonon  energy  above  which  V  =  0, 


Equation  (2-13)  has  a  solution 


dTo) 


(2-14) 


MO)  = 


sinh  1/N(0)V 


For  the  weak  coupling  interaction,  where  N(0)V  <<  1,  the  gap  parameter 
A(0)  is  given  by: 


(2-15) 


A(0)  =  2  exp  (- 

q 


) 


The  BCS  energy  gap  parameter  at  0°K  is  plotted  in  fig.  (2,2c)  as  a  function 
of  energy  e . 

At  higher  temperatures  there  is  always  the  possibility  of  thermal 
excitations  from  the  pair  states.  Because  of  these  excitations,  the 
number  of  possible  pairings  will  be  reduced,  which  reduces  the  binding 
energy  of  the  system  because  some  states  will  not  be  available  for  pairing 
interactions.  This  of  course  reduces  the  gap  parameter  A(T).  In  the  BCS 
theory,  the  energy  gap  is  defined  by  an  implicit  relation  of  the  form 


tahn[(e2  +  A2)2/2  KgT] 

=  /  q  - - - - —  de. 

0 


(ez  +  Az) 


2\  2 


(2-16) 


N(0)  V 
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Fig.  (2-2) 


2 

a  -  Pair  occupation  function  v  as  given  by  equation  (2-9). 

k 

b  -  Quasi-particle  energy  as  given  by  equation  (2-11). 


c  -  The  BCS  energy  gap  A(0)  as  a  function  of  energy  e  . 


c 
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The  transition  temperature  Tc  is  obtained  by  setting  A  equal  to  zero  in 
equation  (2-16) .  In  the  weak  coupling  limit  this  gives 

(2-17)  Tc  =  1.14  eD  exp  (-  )  . 

At  T  =  0°K,  the  half-gap  can  be  written  as 


(2-18)  A (0)  =  1.76  KbTc  . 


The  gap  is  not  expressible  as  a  function  of  temperature  in  an  explicit 
form  over  most  of  the  temperature  range.  Just  below  Tc,  however,  we 
have 


(2-19) 


3.03  (1-t), 


where  t  is  known  as  the  reduced  temperature.  The  gap  has  been  tabulated 
over  the  whole  temperature  range  by  Muhlschlegel  (1959)  from  a  numerical 
calculation.  The  temperature  dependence  of  the  gap  is  shown  in  fig.  (2.3). 


B .  Elementary  Excitations  and  Density  of  States 

The  BCS  formalism  requires  that  the  Bloch  electrons  in  a  metal  be 
replaced  by  quasi-particles,  which  are  excitations  from  the  ground  state. 
These  are  not  exact  eigenstates  of  the  system,  since  the  Hamiltonian  for 
the  system  involves  their  interactions  with  each  other  and  with  other 
excitations,  such  as  phonons  and  plasmons.  However,  they  may  be  treated 
as  approximate  eigenstates  with  a  finite  lifetime. 


Fig.  (2-3)  Temperature  dependence  of  the  energy  gap. 


A/A(oj 
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Bogoliubov  (1958)  developed  an  equivalent  theory  to  BCS.  He  introduced 

another  set  of  operators  formed  from  the  annihilation  and  creation  operators 
* 

Ck*  Ck’  a  canonical  transformation.  This  new  set  is 


Y 


(2-20) 


-> 

k 


* 

It 


u  c 
k  k 


■k 

U  C 

k  -k 


-vc 
k  k 


V  c 
k  k 


Y 


-+ 

k 


Y-S 


u  c 
k  k 


u  c 
k  -k 


+ 


* 

V  c 
k  -k 


* 

V  c 
k  k 


The  vacuum  state  of  these  operators  is  uniquely  defined  by 


(2-21) 


Y  0> 
1 


yJo>  = 
-k 


0 


According  to  (2-20)  y_^  has  the  effect  of  creating  an  electron  in  the 

k 

state  k,  with  amplitude  u  ,  and  at  the  same  time  destroying  an  electron 

k 

-y 

in  state  -k  with  amplitude  v  .  The  probability  amplitudes  u  ,  v  are 

k  k  k 

defined  by  equations  (2-9)  and  (2-10).  If  k  is  well  above  the  Fermi 

2  2 

surface,  so  that  e  >>  A  ,  then  u  =1  and  v  =0.  The  excitation 

k  k  k  k 

will  look  like  an  ordinary  electron.  However,  well  below  the  Fermi 
2  2 

surface,  v  =1  and  u  =0;  the  excitation  now  corresponds  to  the 
k  k  ^ 

destruction  of  an  electron  -  k,  i.e.  to  a  hole.  In  between,  in  the  range 


of  energy  around  A  ,  each  quasi-particle  is  a  mixture  of  an  electron  in  k 

and  a  hole  in  -  k.  The  energy  of  the  quasi-particle  in  the  state  y  |0> 

is  E  ,  and  its  energy  in  the  state  y  y  | 0>  is  2  E  ,  where  E  is 
k  -k  k  k  k 

defined  by  equation  (2-11) .  The  minimum  energy  of  an  allowed  excitation 


is  2  A,  because  excited  single  particles  can  only  be  produced  in  pairs. 


m 
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There  are  no  energy  states  available  in  the  energy  gap.  Therefore 
the  number  of  states  in  an  element  of  energy  dE  is  the  same  as  that  in 
de ,  and  one  can  write 

(2-22)  Ns(E)dE  =  Nn(e)de-; 


where  NS(E)  and  Nn(s)  are  the  density  of  states  of  the  superconducting 
and  normal  phases  respectively.  Using  equation  (2-11),  the  density  of 
superconducting  state  is 


(2-23) 


NS(E)  = 


N  (e) 
n  v 


(E2  -  A2)2  +  A 


dA 

de 


The  BSC  constant  energy  gap  model,  namely, 


A(e)  =  constant  for 
=  0  otherwise, 

simplifies  the  expression  for  Ng(E) 


e  <  hm 

q 


to 


(2-2A) 


N„(E)  =  N(0) 


for  |  E  |  >_  A 


(E2  -  A2) 


2\2 


=  0  for  E  <  A 


where  Nn(e)  is  replaced  by  N(0).  Equation  (2-24)  is  illustrated  in 
fig.  (2.4).  Notice  that  NS(E)  has  infinite  singularities  at  E  =  +  A, 
and  for  E  >>  A,  Ng (E)  =  N(0),  the  density  of  states  of  one  spin  in  the 
normal  metal  at  the  Fermi  energy. 


Fig.  (2-4)  The  BCS  density  of  states  in  a  superconductor. 
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C .  Thermodynamic  Properties 

For  a  metal  in  the  superconducting  state  it  is  customary  to  separate 
the  total  specific  heat,  Cs,  into  a  lattice  contribution  C£s  and  an 
electronic  contribution  C  .  The  lattice  part  is  assumed  to  take  the 
same  values  in  the  normal  and  superconducting  states.  Early  measurements 
of  the  electronic  specific  heat  gave  Ces  «  T3,  in  agreement  with  the  Gorter 
and  Casimir  two  fluid  model.  Deviations  from  the  T3  law  were  first  observed 
by  Brown,  Zemansky  and  Boorse  (1953).  The  BCS  theory  predicts  an  exponenti¬ 
ally  decreasing  electronic  specific  heat  as  T  tends  to  0°K. 

The  entropy  in  the  superconducting  state  may  be  expressed  as 


(2-25) 


[In  (1  +  exp  (-  BE  ) }  +  BE  f  ] 

k  k  k 


* 


where  B  =  : — —  >  and  f  is  the  Fermi  function  for  the  quasi-particles. 

kbt  k 

The  electronic  specific  heat  can  be  expressed  as 

(2-26)  Ces  =  T  §  -  -  6  If-  =  -  4  KBB  [  B  E  -A  . 

k>kF  k 


or 


(2-27) 


'es 


4  Kb62 


l 


->  -> 
k>k. 


->■ 

k 


d-f.) 

k 


[E 


-> 

k 


+ 


£_d_ 
2  dB 


(A2)] 


The  expression  for  C0S  is  interpreted  as  the  specific  heat  due  to  electrons 
and  holes  with  the  modified  spectrum 


E 

-> 

k 


_1_ 

2.  2 


) 


k 


» 


17 


and  the  appropriate  change  in  the  condensation  energy  with  temperature. 
At  sufficiently  low  temperatures  for  which  2A(T)  »  K  Tc,  the  ratio 
Ces/YTc  is  given  by 


(2-28) 


{3 


A(T) 

^  kbtc 


+ 


( Aili\  } 
VKeTcJ 


t 


where  and  x^  are  first  and  third  order  modified  Bessel  functions 

of  the  second  kind.  Equation  (2-28)  can  be  written  in  the  form 


C 

(2-29)  ~r-  -  8.5  exp  (-  1.44  Tc/T) 

-  26  exp  (-  1.62  Tc/T) 


2.5  <  Tc/T  <  6 
7  <  Tc/T  <  11. 


At  the  transition  temperature,  the  jump  in  the  specific  heat,  Ces  -  yT|tc 
is  given  by 

(2-30)  KBN(0)  $2  JTC 


The  derivative  dA2/dg  can  be  obtained  from  the  relation  (2-19),  between 
A  and  T.  The  BCS  calculation  gives  the  value  1.43  yTc  for  the  jump  in 
the  specific  heat,  while  the  Gorter  and  Casimir  model  gives  2  yTc. 

Swihart  (1962)  solved  the  BCS  integral  equation  (2-12)  numerically 
for  non-separable  interactions  of  a  form  more  closely  related  to  that 
derived  by  Frohlich  (1952)  and  by  Bardeen  and  Pines  (1955) .  The  effect 
of  the  Coulomb  interaction  was  also  included.  It  was  possible  to  calculate 
the  specific  heat  and  the  other  thermodynamic  functions  in  the  same  way  as 
in  BCS  theory  with  the  modification  that  the  energy  gap  is  a  function  of 
both  energy  and  temperature.  He  showed  that 
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(2-31) 


(1)  (2) 

^es  _  Ces  ^es 

yTc  -  ytc  YTc 


(1) 

where  Ces  depends  on  the  temperature  dependence  of  the  distribution 

of  quasi-particles  with  constant  energy  gap,  and  goes  to  yTc  as  T 

(2) 

tends  to  T  .  Ceg  depends  on  the  temperature  dependence  of  A(0,T) 
multiplied  by  a  factor  which  is  energy  dependent.  The  factor  was  found 
to  be  one  half  in  the  case  of  weak  coupling  and  the  result  is 


(2-32) 


AC  |  =  _3_  _d_  /a(0,T)\2 

yTc  Tc  2n 2  dt  1  KBtc  J 


In  the  case  of  strong  coupling  the  factor  was  about  0.58,  which  indicates 
that  the  jump  in  specific  heat  at  Tc  is  larger.  However,  this  increase 
is  not  large  enough  to  explain  the  results  for  lead.  In  all  the  treatments 
discussed  so  far,  an  isotropic  energy  gap  model  was  considered.  The 
anisotropy  of  the  energy  gap  will  be  discussed  below. 


D .  Anisotropy  of  the  Energy  Gap 

Anisotropy  of  the  gap  is  well  shown  up  in  ultrasonic  attenuation 
experiments  in  specimens  with  a  long  mean  free  path.  Results  of  experiments 
require  a  half-gap  of  the  form 

(2-33)  A(0)  =  C  KgTc  , 

q  -^g  2  function  of  the  electron  wave  vector  k,  Pokrovskii  (1961) 
proposed  a  model  in  which  he  assumed  an  effective  electron-electron  inter— 
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action  that  depends  only  on  the  direction  of  the  electron  k-vector.  He 
was  able  to  show  in  the  weak  coupling  approximation  that: 

a)  A  change  in  temperature  leads  to  a  change  in  the  magnitude  of 
the  gap  by  a  factor  which  is  independent  of  the  direction. 

b)  The  minimum  value  of  the  gap  at  T  =  0°K  divided  by  the  critical 
temperature  Tc  should  be  less  than  the  BCS  value. 

c)  The  change  in  thermodynamic  quantities  is  qualitatively  the  same 

as  in  the  isotropic  case;  however  the  relative  jump  in  the  specific  heat 

at  T  is  less  than  the  BCS  value, 
c 

The  experiments  show  that  the  relative  jump  in  specific  heat  is  higher 
even  than  the  BCS  value.  Pokrovskii  and  Ryvkin  (1962)  were  able  to  predict 
the  general  effect  of  angular  variations  of  the  energy  gap.  A  more  complete 
theory  was  derived  by  Bennett  (1965)  who  assumed  that  the  strongest  source 
of  gap  anisotropy  is  the  phonon  spectrum.  His  calculation  is  important 
in  suggesting  that  the  temperature  dependence  for  the  gap  agrees  with  BCS. 

The  BCS  expression  for  the  ratio  of  accoustic  attenuation  in  an 
isotropic  superconductor  to  that  in  a  normal  metal  is  given  by 

(2-34)  —  =  2  f  (A/KbT) 

an 

where  f  is  the  Fermi  function.  In  the  anisotropic  superconductor 
Pokrovskii  found  that 

1_ 

(2-35)  —  =  C (T/Tc) 2  exp  (-  Amin/KBT)  , 

an 

where  Amin  is  the  minimum  value  of  the  gap  on  the  curve  defined  by 


- 
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v.q  =  0  on  the  Fermi  surface,  q  is  the  momentum  of  the  sound  wave  and 
v  is  the  Fermi  velocity.  Equation  (2-35)  has  the  same  temperature 
dependence  as  that  of  BCS  theory. 

Another  method  of  observing  the  anisotropy  is  by  means  of  tunneling 
experiments,  which  will  be  discussed  in  the  next  chapter. 


CAHPTER  3 


ELECTRON  TUNNELING 


A.  Tunneling  Theories 


The  most  direct  measurement  of  the  energy  gap  has  been  provided  by 
the  work  with  tunnel  junctions  initiated  by  Giaever  (1961).  A  tunnel 
junction  is  composed  of  two  metallic  films  (1  and  2),  separated  by  a  thin 
insulator.  The  action  of  the  insulator  is  to  produce  a  potential  barrier 
between  the  metallic  films,  wThich  impedes  the  flow  of  electrons  between 
them.  The  barrier  is  assumed  to  be  thin  enough  to  permit  electron  penetra¬ 
tion  by  quantum-mechanical  tunneling. 

The  matrix  element  M  for  a  transiton  of  an  electron  from  a  state 

— ^  ^  * 
kk 

-  y  — ^  ^ 

k  in  metal  1  to  a  state  k  in  metal  2  by  an  applied  voltage  V  is 
assumed  to  be  constant  and  the  same  in  the  superconducting  and  the  normal 
states.  With  this  assumption,  the  probability  of  a  transition  per  unit 
time  is 


(3-D 


k 

1-2 


2n 

IT 


I 

k' 


|M  I2  f (e  )[1 
kk  k 


f (e  +  V)]  6(e  -e  ) 

k"  k  k' 


In  this  expression  f(e  )[1  -  f(e  +  V)]  is  the  probability  that  the 

k  k  + 

y  ~ 

^j^2_^ial  state  k  is  occupied  and  the  final  state  k  is  empty  simultane- 

ously.  The  net  tunneling  current  I  is  then  given  by 


(3-2) 


e  KPl-2  -  P2-l> 
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If  the  density  of  electron  states  is  constant,  and  is  given  by  its  value 
at  the  Fermi  surface  N(0),  when  the  metal  is  in  the  normal  state;  and 
is  equal  to  that  of  the  BCS  expression  (2-24)  when  the  metal  becomes 
superconducting,  then  the  summation  (3—2)  may  be  transformed  to 

oo 

(3-3)  I  =  C  /  dc  p1(e)p2(e+v)  [f(e)  -  f(e+v)]  , 

—  OO 

where 

c  -  f  Imi2I2  Nln(0>  v°>  > 

and  p-^(e)  is  known  as  the  reduced  density  of  states.  It  is  unity  for 
the  normal  metal  and  is  given  by, 


(3-4)  ps 

=  |e  1 

42  +  A2 

for 

|e  |  >>  A 

=  0 

for 

|e|  <  A  , 

for  a  superconductor; 

n  and  s  refer 

to  the 

normal  and  the  superconducting 

states  respectively. 

The  following  three  cases  are  to  be  considered: 


(1)  Inn  ;  both  metals  are  in  their  normal  states. 

The  temperature  of  the  junction  and  the  applied  voltage  are  such  that 
KgT  <<  and  V  <<  Under  these  conditions,  the  tunneling  current  Inn 

can  be  written  as 


(3-5) 


Fig.  (3-la)  illustrates  the  densities  of  electron  states  in  two  normal 
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metals  at  non  zero  temperature  along  with  the  I-V  characteristic. 

(2)  I  ;  metal  1  only  in  its  superconducting  state. 

In  this  case  the  reduced  density  of  states  of  metal  2  is  unity,  and 

that  of  the  superconductor  is  given  by  the  BCS  reduced  density  of  states. 
The  current  I  at  0°K  is  given  by 

0  1 

(3-6)  C  /  pu(e)  de  =  C(V2  -  A^)~  ,  |v|  >_  Aj 

-V 

=  0  ,  |V|  <  Aj  . 

The  limit  of  integration  is  now  from  (-V  to  0) ,  because  of  the  cut  off 
in  the  Fermi  factor  at  0°K,  For  higher  temperatures,  the  solution  for 
the  current  I  corresponding  to  the  potential  region  |v|  <  A  is  given 
by 

00  m-i  /mAj 

(3~7)  Isn  =  2  C  I  t-1)  Xi\kTt 

m=l  o 

where  Xj  is  the  modified  Bessel,  function  of  the  second  kind.  The  density 
of  states  and  the  I-V  characteristic  are  shown  in  fig.  (3.1b).  Notice 
that  for  T  ±  0°K,  there  is  a  small  current  for  V  <  A  due  to  the  thermally 
excited  electrons  which  can  tunnel  into  the  other  metal. 

(3)  I__  ;  both  metals  are  superconductors. 

o  O 

The  most  interesting  experiments  are  those  in  which  the  two  films 
are  different  metals.  A  typical  set  of  investigations  was  that  made  on 

Al-Pb  junctions  by  Giaever  (1961). 

At  zero  temperature,  there  will  be  no  current  until  the  applied 
voltage  corresponds  to  A^  +  A^  .  At  this  voltage,  the  current  increases 
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very  sharply  until  it  approaches  Inn.  When  the  temperature  is  different 

from  zero,  there  will  be  a  net  flow  of  current  from  metal  2  to  metal  1, 

if  a  small  positive  voltage  is  applied  to  metal  1.  As  the  left  edges 

of  the  two  gaps  come  into  coincidence  at  V  =  A  -  A  ,  there  will  be  an 

2  i 

increase  in  current.  This  occurs  because  at  the  temperature  of  the 
experiment  there  are  still  a  considerable  number  of  vacant  states  in 
the  lower  band  of  metal  1  which  would  not  exist  at  0°K.  As  soon  as  the 
voltage  is  increased  further,  the  current  drops  because  there  are  now 
few  states  available.  This  continues  until  a  stage  is  reached  when  the 
right  hand  band  of  1  coincides  with  the  left  hand  band  edge  of  2  at  which 
V  =  A^  +  A^.  Beyond  this,  the  current  starts  to  increase  again.  The 
density  of  states  as  well  as  the  I-V  characteristics  are  shown  in  fig.  (3.1c). 

Cohen,  Falicov,  and  Phillips  (1962)  have  discussed  tunneling  of 
electrons  in  a  superconductor-normal  junction  on  the  basis  of  an  effective 
Hamiltonian, 

(3-8)  H  =  Hn  +  Hs  +  Ht  , 

where  H„  and  Hc  are  exact  Hamiltonians  for  the  normal  and  superconducting 
metals,  respectively,  and  H-p  is  an  operator  which  transfers  electrons  from 
one  metal  to  the  other: 

(3-9)  H  =  l  (T^  C*  +  T  C*  )  . 

-+  ->  kq  ko  ka  kq  qa  ka 

k,q,a 

Here  k  and  q  refer  to  states  in  normal  and  superconducting  metals,  a  is 
the  spin  and  C’s  are  creation  and  annihilation  operators  for  usual  quasi¬ 
particle  states  in  both  metals.  By  making  use  of  the  equation  of  motion, 
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Fig.  (3-1)  The  density  of  states  near  the  Fermi  level  along  with  the 

I-V  characteristics  of 

a)  normal-normal  case 

b)  normal-superconductor  case 

c)  superconductor-superconductor  case. 


T  =  o°K  T  >  oV 


v 


V 
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(3-10)  im  =  [Ns ,  H]  =  [Ng,  Ht]  , 

they  found  an  expression  for  the  time  rate  of  change  of  the  number  of 
electrons,  and  thus  the  tunneling  current  in  the  superconductor. 

When  this  theory  is  extended  to  the  two  superconductor  case,  complica¬ 
tions  arise  because  of  the  omission  of  the  Cooper-pair  transfer  terms 
which  are  present  in  the  microscopic  theory.  Josephson  (1962)  and  others 
have  shown  that  a  proper  treatment  leads  to  additional  terms  in  the  tunneling 
current.  These  terms  give  rise  to  the  so-called  Josephson  effect. 

Schrieffer,  Scalapino  and  Wilkins  (1963)  treated  Up  as  a  small  pertur¬ 
bation  on  the  Hamiltonian 

(3-11)  Hq  =  Hn  +  Hs  . 

They  derived  the  expression 

(3-12)  I  «  /  Np2)(E)  N^,1^  (E-V)de  , 

0 

where 

(3-13)  %(E)  =  N (0)  RP  - - - -  . 

(E2  -  A2)2 

The  dynamic  conductance  dl/dV  of  a  tunnel  junction  is  obtained  by 
taking  the  derivative  of  equation  (3-12) .  The  normalized  conductance 
a  =  (dl/dV) S/ (dI/dV)n  is  found  to  be 

(3-14)  RP  - - - y 

(V2  -  A2 (V) ) 2 


Since  from  equation  (3-13) 
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(3-15)  Nt(V)  =  N (0)  - - - r 

[V2  -  A2(V)]2 

the  measurement  of  o  versus  v  will  give  detailed  information  on  the 
superconducting  tunneling  density  of  states  as  a  function  of  energy.  This 
will  be  discussed  in  the  next  section. 

B .  Phonon  Structure 

BCS  theory  describes  the  superconducting  state  for  electron  energies 
close  to  the  Fermi  surface.  However  in  regions  of  higher  energies  the 
theory  is  not  entirely  in  accord  with  experiment.  This  indicates  that 
the  simplifying  BCS  assumption  of  a  cut  off  for  both  Coulomb  and  phonon 
interaction  at  TiWq  requires  modification.  Bardeen  (1959)  has  pointed 
out  that  the  cut  off  may  be  determined  by  the  lifetime  of  the  quasi¬ 
particles  that  can  be  excited  across  the  gap.  The  effect  of  this  on  the 
pair  interaction  has  been  discussed  by  Eliashberg  (1960),  who  succeeded 
in  extending  Bogoliubov  theory  to  include  both  lifetime  effects  for  the 
excited  electrons  and  an  explicit  phonon  distribution.  Green’s  function 
techniques  have  been  used  to  account  for  lifetime  effects.  It  has  been 
found  that  the  energy  gap  parameter  A  is  a  complex  function  of  energy. 
Swihart  (1962)  and  Morel  and  Anderson  (1962)  assumed  a  model  in  which 
the  interaction  V  is  energy  dependent.  They  showed  that  lifetime  effects 
are  too  small  to  justify  cutting  off  the  Coulomb  repulsion  at  .  There- 

q 

fore  they  include  in  the  interaction  a  repulsive  part  (V  >  0)  for  energies 

si  >  Tim  .  The  resulting  variation  of  the  energy  gap  at  0°K  as  a  function 

k  q 


28 


of  has  been  shown  by  Morel  and  Anderson  to  have  the  form  represented 

in  fig .  (3-2) . 

Equation  (2-24)  indicates  that  any  variation  in  the  energy  gap  should 
be  reflected  in  a  variation  in  the  density  of  states.  Since  the  current 
in  the  tunneling  experiment  is  proportional  to  the  density  of  states,  we 
may  obtain  the  relative  change  in  the  density  of  states  directly  by 
plotting  the  relative  conductance  a  versus  the  applied  voltage  V. 
Giaever,  Hart  and  Megerle  (1962)  found  good  agreement  between  the  data 
and  BCS  theory  for  Sn,  In  and  A£.  However  for  Pb,  definite  divergences 
exist  from  the  BCS  density  of  states.  Adler  and  Rogers  (1963)  observed 
similar  divergences  in  the  case  of  In.  Corresponding  to  these  divergences 
in  the  relative  conductances  versus  energy  curves,  there  are  bumps  in  the 
I-V  characteristics.  It  has  been  suggested  that  these  variations  may  be 
accounted  for  by  using  a  variable  energy  gap  parameter  A(c)  which  changes 
sign  at  e  -  KgO-^. 

Schrieffer,  Scalapino  and  Wilkins  (1963),  have  used  the  tunneling 
formalism  to  investigate  the  structure  in  the  density  of  superconducting 
states  due  to  a  model  phonon  spectrum  for  Pb .  Two  Lorentzian  distributions 
whose  parameters  are  chosen  to  fit  the  neutron  scattering  data  for  Pb  are 
shown  in  fig.  (3-3).  They  also  found  that  the  possibility  of  transferring 
electrons  from  one  film  to  the  other  results  in  a  complex  energy  gap 
parameter,  and  that  the  reduced  density  of  states  is  given  by 


pt(e) 


=  R.P . 


(E2  -  A2)2 


where  A  =  A^  +  i  A^.  Using  a  binomial  expansion  valid  for  A/E  <<  1 
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Fig.  (3-2) 


The  energy  gap  at  0°K  as  a  function  of  energy  as  shown 
by  Morel  and  Anderson  (1962) . 


Energy  gap 
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Fig.  (3- 


The  two  Lorentzian  distributions  whose  parameters  are  chosen 
to  fit  the  neutron  scattering  data  for  Pb . 
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The  variation  with  energy  of  the  ratio  of  density  of 
states  in  superconducting  lead  to  that  in  the  normal  metal.  Curve 
1  is  the  ECS  plot,  curve  2  the  tunneling  measurement.  Also  shown 
is  the  phonon  spectrum  suggested  by  the  experimental  density  - 
of-states  plot. 
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we  find 

PT(E)  =  1  +  -i  • 

So  in  an  energy  region  where  A2  is  larger  than  A1 ,  the  tunneling  density 
of  states  falls  below  that  of  the  BCS  theory. 

Rowell  ej:  ajL  (1963)  looked  for  the  sharp  peaks  in  d2I/dV2  that  occur 
at  voltages  corresponding  to  inflections  in  dl/dV.  These  peaks  are  related 
to  the  Van  Hove  singularities  of  the  phonon  spectrum. 

Van  Hove  (1953)  proved  that  for  a  three  dimensional  crystal,  there 
exist  in  the  phonon  spectrum  (g  (oo))  versus  to  at  least  two  infinite  discon¬ 
tinuities  in  dg/dto,  and  at  the  upper  end  of  the  spectrum  dg/du)  =  -  °°. 

These  infinite  discontinuities  occur  at  the  values  of  to(qc)  for  which 
the  critical  points  q  Ts  are  defined  by  doo/dqL  =  0.  The  effect  of 

C  4c 

the  critical  points  on  the  second  derivative  versus  voltage  plot  of 
N-S  and  S-S  junctions  has  been  considered  by  Scalapino  and  Anderson. 

Their  result  is  that  superimposed  on  the  density  of  states,  is  a  varia¬ 
tion  that  is  attributed  to  the  general  shape  of  the  phonon  spectrum,  with 
singularities  of  the  type  represented  in  table  (3.1). 

Earlier  Morel  and  Anderson  had  pointed  out  that  an  Einstein  phonon 
peak  at  energy  E£  will  produce  structure  in  the  density  of  states  at 
energies 

A  +  nEg  ,  (n  =  1,  2,  .  .  .  .) 

These  types  of  singularities  are  well  discussed  by  Rogers  (1964)  for 


Pb ,  Sn  and  In  samples. 
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Table  (3.1) 


Junction 

dg/dv 

d2I/dv2 

Energy  of  the 
critical  point 

N-S 

00 

logarithmic  singularity 
and  jump  discontinuity 

A  +  Ec 

| 

s-s 

00 

logarithmic  singularities 

2A  +  Ec 

r 

1 

finite. 

CHAPTER  4 


EXPERIMENTAL  TECHNIQUES 

A.  Sample  Preparations 

A  tunnel  junction  consists  of  an  evaporated  base  layer  of  one  metal, 
a  layer  of  oxide  of  this  metal,  and  an  evaporated  metallic  cover  layer. 

The  films  are  deposited  onto  a  glass  substrate.  Tunnel  junction  prepara¬ 
tion  with  a  base  layer  of  Aluminum  is  described  in  detail  by  Adler,  Rogers, 
Khanna  and  Penner.  The  critical  step  in  the  preparation  of  a  specimen  is 
the  oxidation  of  the  base  layer  to  produce  a  barrier  layer  with  specified 
characteristics.  Aluminum  oxide  can  be  formed  even  in  air  at  room  tempera¬ 
ture.  At  elevated  temperatures  and  in  the  presence  of  water  vapour,  the 
rate  of  growth  of  the  Aluminum  oxide  layer  is  much  faster.  It  is  found 
that  the  resistance  of  tunnel  junctions  is  a  function  of  time,  temperature, 
and  the  metal  deposited  for  the  cover  layer.  At  room  temperature,  the 
junction  resistance  increases  very  rapidly  immediately  after  its  preparation 
and  slower  after  that.  At  liquid  air  temperature,  the  resistance  is  almost 
constant  with  time. 

The  junctions  used  in  this  work  are: 

Sn  -  SnO  -  Sn  ,  Sn  -  SnO-  Pb  ,  and  Pb  -  PbO  -  Pb  . 

They  have  different  behaviour  and  require  a  slightly  different  method  of 
preparation.  A  special  cell,  shown  in  fig*  (4-1),  is  used.  It  was 
constructed  for  the  purpose  of  preparing  double  tunnel  junctions  like 
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Fig.  (4-1)  The  cell  which  is  used  for  tunnel  junctions  preparation. 


Sample  Hoider 


Mask 


i 


Rotating  Device 
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Pb  -  PbO  -  Pb  -  PbO  -  Pb .  The  cell  is  provided  with  a  sample  holder 
situated  just  behind  a  rotating  mask  which  is  mounted  in  the  middle  of 
a  furnace.  The  mask  is  a  separate  part,  that  can  be  easily  removed  for 
cleaning  or  replacement.  The  furnace  temperature  ranges  from  20° C  to 
300°C.  It  was  possible  to  control  some  parameters,  like  humidity  and 
oxygen  pressure,  which  might  affect  the  formation  of  the  oxide  layer. 

The  glass  substrates  used  are  pieces  of  clean  microscope  slide  one 
inch  square.  Indium  is  smeared  at  the  positions  of  the  prospective 
electrical  leads.  Usually  the  leads  are  indium  soldered  to  the  films 
after  the  tunnel  junction  is  completed,  but  for  the  tin  junctions,  it 
was  necessary  to  smear  the  indium  solder  on  the  substrate  before  any 
film  deposition  because  the  tin  films  are  destroyed  by  contact  with  molten 
indium.  It  was  found  that  the  tin-indium  contact  often  showed  a  tunnel 
junction  characteristic  at  room  temperature.  This  could  sometimes  be 
destroyed  by  applying  a  high  voltage  across  the  leads  connected  to  the 
film.  The  formation  of  these  tunnel  junctions  is  not  likely  to  occur 
if  the  indium  is  slightly  scraped  to  remove  any  possible  indium  oxide 
just  before  the  film  deposition.  The  subs trate with  indium  smears  is 
shown  in  fig.  (4-2a) .  A  fairly  thick  base  layer  is  evaporated  in 
(10~6  Torr.)  vacuum  from  a  heated  molybdenum  boat  carrying  a  charge 
of  99.999%  pure  tin.  Measurements  of  the  resistances  of  such  films  show 

O 

that  they  are  approximately  500  A  thick.  The  base  layer  is  shown  in 
fig.  (4- 2b). 

Two  methods  of  oxidation  were  tried;  the  glow  discharge  method 
described  by  Miles  and  Smith  (1963),  and  a  thermal  oxidation  method. 


« 
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Barriers  obtained  by  the  first  of  these  methods  were  very  poor.  Metallic 
bridges  were  likely  to  exist  through  the  oxide  layer.  The  current-voltage 
characteristic  of  one  of  these  junctions  displayed  on  the  screen  of  an  oscill¬ 
oscope  is  shown  in  fig.  (4-3).  Notice  that,  for  small  voltages,  the  current 
increases  with  the  applied  voltage,  after  which  many  discontinuities  occur 
in  the  I-V  characteristic.  These  have  been  reported  by  Giaever  (1961)  and 
observed  by  Rogers  (1964) .  The  thermal  oxidation  method  was  used  to  prepare 
all  the  tunnel  junctions  investigated  in  this  work.  The  base  layer  is 
allowed  to  oxidize  for  two  or  three  hours  at  nearly  75°C  in  dry  oxygen  at 
atmospheric  pressure.  The  oxygen  is  dried  by  passing  it  slowly  over  sodium 
hydroxide  and  through  a  coil  of  copper  surrounded  by  a  mixture  of  dry  ice 
and  acetone. 

The  next  step  is  to  evaporate  a  cover  layer  of  tin  or  lead  to  obtain 
the  required  tunnel  junction.  See  fig.  (4-2c) .  The  evaporation  rate  of 
the  cover  layer  was  equal  or  slightly  less  than  that  of  the  base  layer  to 
avoid  barrier  destruction.  The  evaporation  time  was  long  enough  to  produce 
a  continuous  film.  Resistance  measurements  show  that  the  thickness  of 

O 

these  layers  is  approximately  1000  A  . 

About  half  an  hour  after  a  cover  layer  of  tin  was  deposited,  the 
resistance  of  the  junction  was  essentially  zero.  However,  the  resistance 
increases  with  time,  in  contrast  to  Aluminum  junctions,  slowly  at  the 
beginning  and  faster  at  later  times.  As  soon  as  it  reaches  the  required 
resistance,  it  is  mounted  and  immersed  in  liquid  air.  No  further  change 
in  the  junction  resistance  was  observed  during  periods  as  long  as  one  week. 

The  most  stable  junction  is  Pb-PbO-Pb.  These  films  are  also  evaporated 


Fig.  (4-2) 


The  various  steps  of  the  sample  preparation. 
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Fig.  (4-3)  The  I-V  characteristic  of  a  Sn-Sn  junction  which  contains 
metallic  bridges  through  the  oxide  layer. 
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from  a  molybdenum  boat  onto  a  glass  substrate  placed  very  close  behind  a 
fairly  thick  mask  to  avoid  a  shadow  formation.  The  evaporated  films  were 
made  very  narrow  to  decrease  the  possibility  of  metallic  bridges  in  the 
oxide  layer.  Wider  films  result  in  a  tunnel  junction  with  almost  zero 
resistance.  The  base  layer  was  oxidized  in  the  way  described  above,  but 
the  oxide  layer  formation  is  very  sensitive  to  water  vapour.  A  zero 
resistance  oxide  layer  is  always  obtained  if  the  oxygen  contains  a  small 
amount  of  water  vapour. 


B .  Dewar  Description  and  Temperature  Measurements 

The  apparatus  which  is  shown  in  fig.  (4-4)  consists  basically  of  a 
glass  dewar  for  liquid  helium  with  provisions  for  pumping  on  the  helium, 
and  an  outer  dewar  containing  liquid  nitrogen  which  acts  as  a  radiation 
shield  for  the  helium.  A  heater,  consisting  c£  a  500  ft  carbon  resistor 
and  a  small  piece  of  copper,  is  located  at  the  bottom  of  the  internal 
dewar.  Its  function  is  to  remove  the  liquids  from  that  dewar.  The  inner 
dewar  is  provided  with  a  loose  plug  of  styrofoam  at  the  top  to  reduce  the 
evaporation  rate  of  liquid  helium. 

Samples  are  immersed  directly  in  the  helium  bath;  this  provides  a 
good  thermal  coupling  between  the  bath  and  the  sample.  Temperatures  are 
determined  from  measurements  of  the  helium  vapour  pressure.  The  pressure 
is  controlled  by  means  of  a  manostat  with  a  flexible  membrane.  In  such 
an  arrangement,  temperature  inhomogeneities  in  the  liquid  helium  may  make 
the  recorded  vapour  pressure  a  rather  poor  indication  of  the  temperature 
at  points  within  the  liquid.  A  stationary  temperature  can  be  maintained 
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Fig.  (4-4)  The  glass  dewar  for  liquid  helium. 
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if  the  corresponding  pressure  is  approached  by  reducing  the  pressure 
slowly  and  regularly.  Small  corrections  to  the  manometer  readings  for 
temperature  and  g-value  differences  from  standard  values  have  been 
neglected.  The  accuracy  in  temperature  was  about  +0.1  per  cent  in 
the  temperature  range  from  2  to  4°K. 


C .  Instrumentation 

The  I-V  tracing  circuit  is  shown  in  fig.  (4-5).  It  has  been  designed 
and  constructed  by  Rogers  0.967).  The  circuit  is  provided  by  a  sensitive 

current-detector  and  a  small  impedence  current  source  (nearly  vertical 

load  line),  in  order  to  detect  currents  as  small  as  10”7  amperes.  The 

idea  of  the  circuit  is  shown  in  fig.  (4-6).  The  point  of  intersection  V 
satisfies  the  equations: 

VB  ” * * * *  V 

I  =  - - —  for  the  load  line 

and  I  =  I(V)  for  the  tunneling  current. 

The  locus  of  the  point  V  as  the  load  line  moves  outward  with  a  constant 
slope  yields  the  I-V  characteristics  of  the  tunnel  junction.  The  supplied 
voltage  Vg  can  be  swept  smoothly  and  uniformly  in  time  at  a  rate  which  is 
polarity,  range  and  vernier  adjustable.  The  sweep  amplitude  control 
attenuates  the  output  from  the  sweep  generator,  shown  in  fig.  (4-7),  to 
a  convenient  range  for  energy  gap  measurements. 

The  a  versus  v  traces  are  obtained  by  means  of  a  conductance  bridge 
described  by  Adler,  Roger  and  Woods  (1964).  The  bridge  has  been  modified 
and  is  described  by  Rogers  (1967) . 


_ 


Fig.  (4-5)  The  I-V  tracing  circuit. 
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Fig.  (4-6)  Schematic  showing  the  idea  of  the  I-V  tracing  circuit. 
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The  sweep  generating  circuit. 
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CHAPTER  5 


EXPERIMENTAL  RESULTS  AND  DISCUSSION 
A.  The  Relative  Jump  in  Specific  Heat 

Fig.  (5-1)  shows  a  typical  I-V  characteristic  for  a  Sn-Pb  tunnel 
junction  at  3°K  along  with  a  sample  calculation  of  the  energy  gaps. 

(The  break  in  the  characteristic  within  the  negative  resistance  region 
is  attributed  to  the  onset  of  oscillation) .  Two  sets  of  the  I-V  charac¬ 
teristics,  at  different  temperatures,  of  one  of  the  Sn-Pb  junctions  and 
one  of  the  Sn-Sn  junctions  are  shown  in  figs.  (5-2)  and  (5-3).  Notice 
that  the  Sn  films  still  show  superconductivity  even  at  3.79°K,  while 
the  transition  temperature  of  bulk  tin  is  3.73°K.  A  value  3.8°K  was 
chosen  for  Tc  and  used  to  plot  t^sn(T) /KTC] 2  versus  t  =  T/Tq  in  the 
range  0.9  <  t  <  1.  These  plots  are  the  straight  lines  shown  in  figs. 

(5-4)  and  (5-5).  The  actual  values  of  Tc  were  determined  by  adjusting 
the  slope  of  the  lines  to  intersect  the  t  -  axis  at  t  =  1.  The  transition 
temperatures  were  3.81°K  for  Sn  in  a  Sn— Pb  junction  and  3.83  K,  3.84  K 
for  Sn-Sn  junctions.  At  these  temperatures,  a  value  1.61  +  0.02  for  the 
average  of  the  relative  jump  m  specific  heat  was  obtained  from  the  slopes* 

of  the  lines  and  equation  (2-32). 

Tunnel  junctions  of  Sn— Pb  with  resistances  from  0.3  to  4  ohms  at 

O 

liquid  nitrogen  temperature  and  thickness  in  the  range  from  420  A  to 
3670  A°  were  investigated  by  Shigi  (1968).  The  transition  temperature 
for  Sn  in  all  the  samples  was  found  to  be  3.86°K,  in  good  agreement  with 
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Fig.  (5- 


1)  The  typical  I-V  characteristic  of  a  Sn-Pb  junction  at 


3°K  along  with  a  sample  calculation 


of  the  energy  gaps. 
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Fig.  (5-2)  The  I-V  characteristics  of  a  Sn-Pb  junction  at  different 


temperatures . 


Sn-Pb  Junction 
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Fig.  (5- 


The  I-V  characteristic  at  different  temperatures  of  a  Sn-Sn 
junction. 
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Fig.  (5-4)  (A/kTc)2  versus  t  =  T/Tc  for  a  Sn-Sn  junction. 
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Fig,  (5-5)  (A/KTC)2  versus  t  =  T/Tc  for  Sn-Pb,  and  Sn-Sn  tunnel  junctions. 
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the  values  determined  in  this  work. 

Lock  (1951)  has  associated  the  increase  in  Tc  of  tin  films  above 
the  value  for  the  bulk  tin  with  stress  in  films  constrained  to  follow 
the  thermal  contraction  of  a  glass  substrate.  Jennings  and  Swenson  (1958) 
showed  that  a  compressive  stress  decreases  the  critical  temperature  while 
tension  increases  it.  Tin  films  deposited  onto  glass  substrates  were  also 
investigated  by  Blumberg  and  Seraphin  (1962) .  From  resistance  measurements 
they  found  that: 

a)  The  transition  temperature  varied  considerably  from  sample  to  sample, 
even  when  the  thickness  and  all  known  control  parameters  are  reproduced. 

b)  There  is  a  regular  increase  in  Tc  with  decreasing  film  thickness, 

O 

for  film  thickness  less  than  300  A  .  c)  Much  of  the  variation  in  Tc 
at  the  same  film  thickness  is  contributed  by  anisotropy  effects.  The 
crystalline  orientation  of  the  film  is  significant  in  determining  the 
magnitude  of  the  stress  and  the  variation  in  critical  temperature  relative 
to  bulk  tin.  d)  The  transition  width  increases  with  decreasing  thickness 

o 

until  for  70  A  ,  the  width  is  of  the  order  of  0.070  K.  This  is  in  contrast 
with  the  result  of  Walmsley  and  Campbell  (1967),  that  the  current  step  in 
the  current— voltage  characteristics  of  a  tunnel  junction  is  broader  in 
thick  films  than  in  thinner  ones.  In  order  to  minimize  the  effects  of 
thermal  stress,  the  films  should  be  deposited  on  substrates  whose  thermal 
expansion  coefficients  closely  match  those  of  the  film  material. 

Accurate  determination  of  the  electronic  specii_ic  heat  Ces  can  be 

*  Rogers  claimed  that  the  transition  temperature  determined  from 
the  resistance  of  the  film  itself  is  unreliable  due  to  edge  effect.  In 
Blumberg ’s  work,  the  films  are  ruled  to  remove  the  edge  effects. 
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made  from  calorimetric  measurements  only  on  those  superconductors  for  which 
Ces  is  a  large  fraction  of  the  total  specific  heat.  These  measurements 
provide  a  test  of  the  assumption  that  the  lattice  contribution  to  the 
specific  heat  is  the  same  in  both  superconducting  and  normal  phases. 

At  very  low  temperatures,  this  is  in  conflict  with  the  experimental 
results  for  some  superconductors.  Brayant  and  Keesom  (1960)  found  that 
Cs  is  less  than  C^n  below  0.8°K  for  Indium.  As  a  consequence  of  this 
behaviour,  they  concluded  that  it  is  not  possible  to  estimate  Ces  directly 
from  specific  heat  data.  The  electronic  specific  heat  as  a  function  of 
temperature  at  very  low  temperatures  can  not  be  obtained  from  energy  gap 
measurements,  since  the  gap  is  almost  constant  over  the  range  0  <  t  <  0.4. 
However,  close  to  the  transition  temperature,  the  energy  gap  is  strongly 
temperature  dependent.  This  provides  a  measure  of  the  jump  in  Ceg  at  Tc. 

This  is  the  first  determination  of  the  relative  jimp  in  Ceg  for  Sn 
by  measuring  the  energy  gap  from  the  I-V  tunneling  characteristics.  The 
results  are  in  good  agreement  with  the  value  1.6  obtained  from  calorimetric 
measurements.  The  magnitude  of  this  jump  depends  upon  the  strength  of  the 
attractive  interaction,  the  phonon  spectrum,  and  the  band  structure  of  the 
superconductor.  In  Pb  which  is  strong  coupling  and  which  has  a  less 
anisotropic  Fermi  surface  and  phonon  spectrum  than  Sn,  the  magnitude  of 
the  jump  is  2.65  from  calorimetric  measurements  and  2.15  from  tunneling 
measurements  (Frank  and  Keeler  1968) . 


B.  The  Energy  Gap  2A(T)  and  the  Ratio  2A(0)/KTc 


The  ratio  A(T)/A(0)  as  a 


function  of  the  reduced  temperature  t  along 


. 

■J, 
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with  the  BCS  curve  is  shown  in  fig.  (5-6).  For  this  figure,  the  value  of 
2A(0)  has  been  adjusted  to  give  a  good  visual  fit  with  the  BCS  curve 
represented  by  the  solid  line.  From  this  chosen  value  of  2A(0)  -  1.22 
mv,  one  obtains  a  ratio  3.68  for  the  energy  gap  2A(0)/KTC.  The  gap  2A(T) 
also  extrapolates  to  1.22  mv  at  0°K.  In  table  (5-1)  the  values  of  Tc 
and  2A(0)/KTC  determined  in  this  work  are  shown  together  with  some  other 
results . 

The  temperature  dependence  of  the  energy  gap,  from  ultrasonic  attenua¬ 
tion  experiments,  is  shown  in  fig.  (5-7)  after  Morse  (1957).  In  fig.  (5-8) 
we  plot  the  deviations  of  the  results  of  the  electron  tunneling  experiment 
and  the  ultrasonic  attenuation  experiment,  from  that  of  the  BCS  theory 
in  terms  of  the  reduced  temperature  t.  The  average  deviation  <dt>  of 
the  ultrasonic  attenuation  data  from  the  BCS  theory  is  -  5%  for  °*3  <  tBCS 
<  0.7  and  +  5%  for  0 . 7  <_  tgCg  <  1,  while  the  electron  tunneling  points 
deviation  is  -  1%  for  0.4  <_  tg^g  <  0.6  and  +  1.5  %  for  0.6  <_  tg^g  <  1. 

As  mentioned  in  chapter  (2) ,  the  ultrasonic  attenuation  experiment 
requires  A(0)  =  BKTQ  instead  of  the  BCS  expression,  where  B  is  a  function 
of  the  crystallographic  direction.  The  values  of  B  in  the  case  of  Sn 
for  different  orientations  are  given  in  table  (5.2)  below,  after  Shepelev 
and  Filimonov  (1965).  Privoroboskii  (1963)  showed  in  a  theoretical  analysis 
that  the  minimum  energy  gap  of  superconducting  tin  has  not  been  determined, 
and  that  it  is  not  located  along  one  of  the  principal  crystallographic 
directions.  The  smallest  value  2.7  for  B  has  been  determined  by  Zavaritskii 
(1963)  from  electron  tunneling  measurements,  by  using  a  single  crystal  of  Sn 
on  one  side  of  a  tunnel  junction.  Thin  evaporated  films  are  polycrys tallme. 


. 


Fig.  (5-6) 


The  temperature  dependence  of  the  energy  gap. 


0  8 
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Fig,  (5-7)  The  temperature  dependence  of  the  energy  gap  obtained  by 
Morse  (1957)  from  ultrasonic  attenuation  data. 


A(T)/a  (0) 


Fig.  5.7*  The  ratio  J(T)/J( 0)  versus  TITe.  The  solid  curve  is  the  theoretical 
result  of  BCS.  The  points  are  the  data  of  Morse  and  Bohrn  (M21)  from  the  atte¬ 
nuation  of  longitudinal  sound  waves  in  tin. 
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Fig.  (5- 


8)  The  deviations  of  the  results  of  the  electron  tunneling 
experiment  and  the  ultrasonic  attenuation  experiment 
from  that  of  the  BCS  theory  in  terms  of  the  reduced 


temperature  t. 


3C9 


Table  (5.1) 


Tc 

2A(0)/KTC 

Reference 

3.84 

3.68 

This  work 

3.812 

3.65 

Douglass  and  Meservey  1964 

3.73 

3.61 

Rogers  1964 

3.8 

3.51 

Townsend  and  Sutton  1962 

3.85 

3.67 

Walmsley  and  Campbell  1967 

3.86 

3.57 

Shigi  1968 

Table  (5.2) 


Orientation  of 
accous tic 
vector  k 

Maximum 
frequency  MC 

2A0/KTC  =  B 

k  |  (101) 

257 

3.9 

K  |  (111) 

252 

4.8 

K  |  (112) 

257 

4.4 

K  |  (211) 

276 

3.9 

K  |  (301) 

257 

4.1 

K  |  (113) 

257 

4.0 

K  |  (311) 

215 

4.3 
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For  tunnel  junctions  composed  of  such  films,  it  is  impossible  to  observe 
the  anisotropy  in  the  energy  gap  directly.  Although  the  anisotropy  is 
not  completely  removed,  what  one  can  observe  is  the  average  of  anisotropic 
energy  gaps  over  the  area  of  the  tunnel  junction.  Such  an  average  may  be 
approximated  by  assuming  that  all  crystallographic  directions  contribute 
equally.  The  discrepancy  in  the  ratio  A(0)/KTc  of  film  tunnel  junctions 
depends  upon  which  crystalline  orientation  lies  in  the  plane  of  the  film 
and  which  one  is  perpendicular  to  it.  The  crystallographic  directions 
are  not  reproducible  for  different  films  evaporated  under  conditions 
which  outwardly  appear  to  be  the  same. 


C.  Phonon  Structure  in  Sn 

The  dl/dV  versus  v  characteristic,  at  a  small  sensitivity,  for  a  Sn- 
Sn  junction  at  approximately  2°K  is  shown  in  fig.  (5-9).  In  fig.  (5-10), 
the  same  characteristic  is  shown  at  higher  sensitivity;  the  solid  line 
represents  the  normalized  curve,  whereas  the  dashed  ones  are  not  normalized. 
We  located  the  points  at  which  the  second  derivative  d2I/dV2  would  have 
maxima  and  minima.  These  are  represented  by  the  arrows  in  fig.  (5-10). 

The  maxima  occur  at  the  critical  points  3.3  +  0.1  mv  and  17.6  +  0.1  mv. 

The  minima  occur  at  4.5  +  0.1  mv  and  16.1  +  0.1  mv.  There  may  also  be 
weak  structures  at  7 . 9  +  0 . 1  mv  and  10 . 8  +  0 . 1  mv.  The  different  types 
of  singularities  are  not  analysed  here,  since  the  d2I/dV2  versus  v  traces 

have  not  been  obtained. 

The  structure  in  tin  is  complicated  and  it  was  hard  to  locate  the 


values  of  energy  at  which  the  structure  occurs 


from  the  first  derivative 


The  dl/dV  versus  v  characteristic  of  a  Sn-Sn  junction  at 

2°K, 


Fig.  (5-9) 


(arbitrary  units) 


— I  > 

ia|  T3 


Sn-Sn  (2°  K) 


1 _ __i 


-2-1  0  1  2 

v(mv) 
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Fig.  (5-10)  The  locations  of  the  critical  points  at  which  the  phonon 
structure  occurs  in  the  dl/dV  versus  v  characteristic. 


Sn-Sn  (2°K) 


(SiiNn  AHvaiiayv) 


61 


versus  energy  plot.  Clearer  analysis  can  be  obtained  from  the  second 
derivative  versus  energy  curves.  This  has  been  done  by  Rowell  and 
Kopf  (1965)  on  Sn-Sn  junctions.  They  found  that  the  tin  phonon  spectrum 
extends  to  17.7  mv,  with  critical  points  as  low  as  3.4  mv.  Rogers  (1964) 
observations  using  Al-Sn  junctions,  are  in  good  agreement  with  the  work 
of  Rowell  and  Kopf.  He  observed  further  that  the  tin  spectrum  has  a 
broad  transverse  group  at  about  4.5  mv,  and  a  longitudinal  group  at  16 
mv.  The  results  obtained  from  this  work  are  in  good  agreement  with  Rowell 
and  Kopf,  and  Roger’s  results. 


D.  Phonon  Structure  in  Pb 

The  phonon  structure  in  Pb  is  simpler  than  that  of  Sn.  In  a  Pb-Pb 
junction,  the  energy  gap  at  4.2°K  is  found  to  be  2.56  mv,  obtained  from 
the  a  versus  v  trace  shown  in  fig.  (5-11).  In  fig.  (5-12a)  ,  we  observe 
two  transverse  groups ;  the  first  is  located  around  3.7  +  0.1  mv,  and  the 
second  at  4.5  +  0.1  mv.  A  longitudinal  group  is  observed  at  8.5  +  0.1  mv, 
while  some  other  structures  are  found  at  1.7  +  0.1  mv,  3.1  +  0.1  mv  and  9.0  + 
0.1  mv.  The  signs  and  the  location  of  the  critical  points  are  in  a  fairly 
good  agreement  with  Rowell  and  Kopf  results  (1965)  and  with  Brockhouse  et  al 
(1962)  and  Rogers  (1964).  Rowell  and  Kopf  observed  structures  at  1.75,  3.05 
and  9.0  mv,  a  transverse  group  at  4.4  mv  and  a  longitudinal  group  at  8.5  mv. 
The  structure  at  9.0  mv  was  also  observed  by  Brockhouse  as  a  longitudinal 
group  along  different  crystal  directions.  We  believe  that  the  structure  at 
9.0  mv  is  the  second  harmonic  of  that  at  4.5  mv,  and  the  one  observed  at  17 
mv  in  fig.  (5-12b)  is  also  the  second  harmonic  of  that  at  8.5  mv. 
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Fig.  (5-11)  The  energy  gap  in  Pb  calculated  from  the  dl/dV  versus  v 

characteristics  of  a  Pb-Pb  junction. 
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CHAPTER  6 


CONCLUSION 

The  tunneling  techniques  is  a  powerful  method  of  investigating  the 
relative  jump  in  the  specific  heat  at  Tc,  especially  in  the  case  of 
metals  for  which  Ces  is  a  small  fraction  of  Cs. 

Although  the  theories  of  the  energy  gap  anisotropy  predict  that 
the  relative  jump  in  the  specific  heat  should  be  less  than  that  predicted 
by  the  BCS  theory;  the  experimental  result  is  not  in  agreement  with  the 
theories  in  this  respect. 

The  discrepancies  in  the  ratio  2A(0)/KTC  for  the  same  superconductor 
found  from  different  electron  tunneling  experiments  is  almost  certainly 
due  to  thermal  stresses  and  anisotropy  effects. 

The  phonon  structure  results  for  Sn  and  Pb  are  in  good  agreement 
with  results  published  in  the  literature. 


APPENDIX 


TUNNEL  JUNCTIONS  AS  THERMOMETERS 

The  current-voltage  characteristics  possess  direct  information 
about  the  temperature  of  the  metallic  films,  in  their  superconducting 
states,  of  a  tunnel  junction. 

At  a  temperature  T  <  TC°K,  the  thermal  excited  electrons  and  holes 
can  tunnel  from  one  metal  to  the  other  when  a  voltage  V  <  2A  is  applied. 
Rogers  (1964)  has  calculated  the  normalized  current  I  =  Iss/Inn  for  an 
applied  voltage,  V  =  2A;  slightly  less  than  the  energy  gap  2A.  He  found 
that,  for  a  symmetrical  tunnel  junction,  the  relation  between  the  norma¬ 
lized  current  and  the  temperature  can  be  written  as 

1_ 

_  -A  2  -A/T 

A-l  I(V  =  2A)  =  Igg/inn  =  1 *33  (f)  e 

Although  he  introduced  an  approximation,  for  the  BCS  function  "GM, 
in  the  expression  of  the  tunnel  current  integral,  the  experimental  results 
of  our  work  are  in  good  agreement  with  his  calculation. 

The  observed  values  of  Igs/^-nn  are  plotted,  in  fig.  (A— 1)  for  two 
samples  of  Sn-Sn  as  a  function  of  A^T--  along  with  equation  (A-l).  In 
the  same  figure,  other  results  have  been  obtained  from  one  of  Giaever's 

publications  (Giaever ,  et  al  1961)  by  Rogers  (1964). 

The  discrepancy  in  our  results  may  be  due  to  small  difference  in 
energy  gaps  of  the  two  tin  films  of  the  tunnel  junction. 
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Fig.  (A-l) 


The  normalized  current  ISs/^-nn  versus 


A 

T 


along  with  the 


mathematical  expression  (A-l) „ 
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